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For high-frequency scattering by a plane electromagnetic wave incident upon a convex surface, 
the current distribution in the transition region between the illuminated and shadow region, and beyond, 
was considered by V. A. Fock. The effect of the surface being flat and no longer parabolic in the 
penumbra region is considered here. There are two cases considered, one for which a plane wave 
is incident along a plane perfectly conducting surface towards a drop in this surface, this drop rep- 
resented by a locally convex parabolic cylinder smoothly joined to the plane section. An expression 
is derived for the fields on the convex surface in the penumbra region, and deep in the shadow region 
the appropriate creeping wave expression is obtained. For the other case, a plane wave is incident 
towards a convex parabolic section smoothly joined to a planar section with the parabolic section 
lying in the illuminated region, and the direction of the incident wave being along the planar section. 
It is shown that far along this flat section, the total field on the surface is comprised of the incident 
field plus a traveling wave field which decreases as the square root of the distance from the illuminated 
side. 

1. Introduction 

For high-frequency scattering by a plane electromagnetic wave incident upon a smooth convex 
(locally parabolic) surface, the current distribution in the transition region between the illuminated 
and shadow regions was obtained by V. A. Fock [1946]. Fock gave, as an estimate of the width 
of the penumbra region, the following relation 

where R is the radius of curvature at the geometric optics shadow boundary. The current dis- 
tribution was expressed in terms of a contour integral involving Airy integrals, and which, asymp- 
totically, gave the physical optics current distribution in the illuminated region and a creeping 
wave type of current distribution in the shadow region. 

The purpose of this paper is to extend the results to the case where the geometry of the surface 
is cylindrical and is such that over a portion large compared to wavelength, the surface is approx- 
imated by a flat section smoothly joined to a parabolic section. Hence the surface, besides pos- 
sessing a region of infinite radius of curvature, also possesses a discontinuity in curvature, i.e., 
a discontinuity in the second and all higher derivatives of the surface profile. The incident field 
is assumed to be a plane wave traveling along in a direction parallel to the plane of the flat portion, 
and thus the shadow boundary and join coincide. 

There are two cases treated. First considered is the case where the plane wave is propagating 
along the flat portion towards and normal to the junction of the flat section with the parabolic sec- 
tion. The convex parabolic section is thus in the geometric optics shadow zone. Next considered 
is the case where the incident wave is traveling in the opposite direction such that the parabolic 
section is in the geometric optics illuminated region. 

The analysis is restricted to the important case where the polarization of the electromagnetic 
field is perpendicular to the surface. Hence since, in addition, the geometry is cylindrical, the 
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problem reduces to a scalar problem. A briet discussion taking into account the other polariza- 
tion and a more general noncylindrical surface is given at the end. Briefly, the analysis will pro- 
ceed as follows: The unknown field distribution on the surface will be expressed in terms of an 
integral equation which is then simplified by the following process. Physical postulates are placed 
upon the behavior of the unknown field distribution in the limit of small wavelengths. A high- 
frequency asymptotic analysis then reduces the integral equation to a simplified form which is 
amenable to solution with employment of the Laplace transform. The field distribution is easily 
obtained and in the first case considered, is expressed in the form of a contour integral involving 
Airy integrals. Deep in the shadow region this integral formulation is then expressible in terms 
of a sum of residues, physically representing the phenomenon of creeping waves. In the second 
case considered, it is shown that on the flat section, a long distance from the junction, the field is 
comprised of the incident wave plus a traveling wave whose magnitude decreases as the square 
root of the distance of the field point from the join. 

2. Convex Section in the Shadow Region 

A cartesian coordinate system will be used and oriented such that the z-axis lies along the 
junction of the plan and parabolic sections, and the plane x = 0, y< describes the flat portion 
of the surface. The convex section will be in the half-space y > 0, and will be represented locally 
in the neighborhood of the junction by the parabolic cylinder 

*=-^f a) 

where R is the radius of curvature of the parabolic section at the junction. Since the analysis 
will be concerned with the neighborhood of the junction, the curvature of the convex surface will 
not be prescribed for large y. 

A plane wave will be considered incident along the plane section towards, and normal to the 
junction {x = 0, y=0). The wave will be polarized perpendicular to the conducting surface, with 
the magnetic intensity in the z direction. The incident field has the prescribed form 

H = \ z e ik v 



E=-i x e ik yVJjLoleo 



where harmonic time dependence e~ io)t is assumed. The problem reduces to that of a scalar prob- 
lem involving the magnetic field component, which must satisfy the following integral equation 
on the surface [Maue, 1949]: 

dn 2tt J driQ 

with 



I = 0, iiOP) = 2u (P) + ± j |£ U(Q)ds Q ( 2 ) 



e ikp 

G = and p= |xp — xq|. 

Since the source of illumination is on the surface, Maue's equation implies that 

2uo = e ik v. 

Using the results of (9), (10), and (11) of [Weston, 1962], and therein setting a x = 0, and a 2 =l/R, 
(2) reduces to 



u{y) = e iky- l A f 30 dt^m\kr)w{t, y) 
£ J— oo r 



(3) 
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where 



(y-t) 2 + 



±{Zt 2 -af) 2 \ 



w(t, y) = at(t — y) — - (at 2 — ay 2 ) 



(4) 
(5) 



with the relationships 



£<0, 

y<0, 



a = 0, 
a = 0, 



J>0, 

y>o, 



a = R~ l ; 
a = RK 



To simplify the analysis and to reduce the integral equation to a form amenable for solution 
the following physical postulates will be placed upon the behavior of the unknown surface field 
u(y). On the plane section the total field will be comprised of the incident wave plus a wave 
reflected by the junction, as indicated by the following relationship: 



y < 0, u(y) = e ik v + e^I-iy). 



(6) 



In addition in the high-frequency limit, it will be assumed that the factor I-(y) is a slowly varying 
function compared to the oscillatory factor exp {—iky}. On the parabolic section the field is a 
wave traveling away from the junction, having the form 



y > 0, u(y) = e iks I+{y), 



ith 



5 = 



Jo 



+ y 2 IR 2 ] 1 ' 2 dy 



(7) 



(8) 



being the arc-length along the surface, measured from the junction. Again it will be assumed 
that the factor I+(y) is a slowly varying function in comparison with the oscillatory factor. In the 
limit where kR^&, the following must hold: /_(y)— >0 and /+(y)— »1. 

From (3), (4), and (5) the following integral equation holding for y < is obtained 



e ik H-(y) 



_\k_ f- 
4/? Jo 



dtI + (t)e mt) HP(hr) 



(t 2 -2ty) 



(9) 



with 



r 2 =(y-t) 2 +^i 



for y > 

e mv)l + (y) = e av-i^ f° —[e ik ' + e~ ikl I-(t)]H[ l \kr) 
4a J-os r 



4« 2 



(9a) 



ik 



[* ^e^I+itm'Kkr) (t - y) 2 (10) 

Jo r 
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^here 



^=(y-^) 2 + ^i (ii) 



and 



(y 2 — t 2 ) 2 



in the first and second integrals respectively. 

Since I+(y) is a relatively slowly varying function, (9) can be written in the form 

J-(y) = 2^R J exp [**/M'» y) dkt 

where v(t, y) is a slowly varying function and 

f(t, y) = y + s{t) + r(y, t). 

Since kj{t 9 y) has no stationary point in the region of integration, the dominant part of the integral 
arises from the lower limit. Thus on integrating by parts one obtains 



'-(r)~o(^)/ + 



(0). 

But on taking the limit as y—> + in (10) one obtains 

/+(0)~1; 

thus I-(y) can be neglected, which is to be expected since the surfaces are smoothly joined. 

For the region y > 0, that portion of the integral in (10) that contains the reflected wave has no 
stationary point; hence its dominant part arises from the upper limit of integration. Since in 
addition the reflected wave is the order of (1/kR) of the incident wave, this integral is small com- 
pared to that portion of the integral in (10) containing the incident field. The dominant portion of 
the last integral in (10) arises from the region ^ t ^ y. The reason is that there exists a stationary 
point in the region ^ t ^ y, and no stationary point in the remaining interval y < t < °°. Hence 
in the high-frequency limit, the integral over the range y < t < o° may be evaluated by integration 
by parts, the dominant portion arising from the lower limit of integration. Thus this integral is 
the order of (kR)' 1 I+(y) and can be neglected compared to the left-hand side of (10). 

With these considerations (10) can be approximated by 

7 lev 2 f ° p ikt ik fy dt 

e ik. I+br) = etk y-i^.j ^£_#(i W _i|^ °L e im) I4t)H (» {kr)(t - y) 2 (13) 

where the r's are specified accordingly by (11) and (12). 

To obtain the appropriate asymptotic solution to (13) for large kR we shall set 

m = (kRyl s (14) 

ky= m 2 £ (15) 

kt=m 2 T (16) 

hiy) = M)- (17) 
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Multiply both sides of (13) by the factor exp {— iks{y)}. In obtaining the appropriate asymptotic 
value of the resulting integral corresponding to the first integral on the right-hand side of (13) the 
following asymptotic expressions are used: 

%-*(y)]~4^ + °(i) (18) 

kr-it-rW + ^m-T) (19) 

where r is specified by (11). Hence in the limit m^> °°, the following integral 

;j,2v f° rlt 

J-» r 



4R 



reduces to 



f • expt{gV[8(g-T)]} J^ l4 f * 

J- tf-0 3 ' 2 V^ Jo 



ith 



x-gr- 



In evaluating the asymptotic value of the remaining integral, the following relations are used: 

fcr ~ (f — T)m? H 

k[s(t) - s(y) + r] ^=j^ + (^ 

where r is specified by (12). Thus in the limit m — > «, the second integral given by 



— iA: 
4fl 



f " - h(t)H\'\kr){t-yf exp {««*)- *(y)]} 
Jo r 



reduces to 



4 V<77 



| |* dr (f - r) 1 ' 2 /(t) exp [- H{ - r) 3 /24] + (^j • ( 2i) 



Using these results (13) reduces to the form 



9 f* ,,-177-/4 /O 

y(f)~exp(-if 3 /6)[l-- 7 =e- i ' r/4 e iz2 dz\-^—- 
V^ Jo 4 



|J dr {£ - r) 1 ' 2 7(r) exp [ - i({ - r) 3 /24] . (22) 

This resulting equation is of a form which is amenable to solution by well known methods. As a 
preliminary the following Laplace transforms will be defined, 

J(p)*°rer*Mde, (23) 
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M (p) = f X e-* e 12 exp ( - ^ 3 /24)<#, (24) 

Jo 

N(p)= [ X exp(-p£-^ 3 /6)[l 2 —e~^ P e^dz]^- (25) 

Jo V7T Jo 

With these definitions, the solution to (22) is obtained by taking the Laplace transforms of 
both sides, yielding 

J(p) = N(p)-- yj-M(p)J(p)- 26) 

Rearranging the terms gives the resulting expression for J (p), 

f p-itt/4 /o~ 1-1 

/(p) = #(p)/|l+§ V^ M(p) J " < 2? ) 



Fortunately M(p) has already been evaluated [Weston, 1962] and has the form 

pi37r/4 



i 



M(p) = 2tt V7t8 \Ai (A[ - i B{) + g 



(28) 



where At and Z?i are Airy integrals [Miller, 1946], with argument — ip2 1/3 . The denominator of 
(27) may be simplified by using the Wronskian relationship of the Airy integrals, becoming 

27riAM[-iB[). (29) 

Unfortunately Nip) has to be computed; for purposes of evaluation, the limits of integration over 
the z variable in (25) will be changed from and x to x and °° , where in doing so use is made of 
the relationship 



J>*-^ 



N(p) reduces to the form 



with 



e i " 2 dz= x Hze"l 4 ; 



N(p)= =_e-W F(p), (30) 



F(p)= \ . e xp(-pi-if/6) [" exp (iz 2 ) ofe <#. 
Jo J X 



Expression F(p) is evaluated explicitly in the appendix. 

An explicit expression can now be obtained for J(p). From (27), (29), (30), and (A. 21) (in 
the appendix) the following expression is obtained: 

J(p)=- i2 4 ' 3 [Va + y 2 J" (Ai(x) - i Bi(x)] dx] 

Ai(q)-iB' t (q) < 31 > 
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where 

q = -ip2^. (32) 

Now J(p) can be obtained by taking the inverse transform yielding 

M) = ~-.\ C+ "e(r>J(p)dp- (33 ) 

Changing the variable of integration by setting p = it2~ 113 and using the Airy integral notation 
used by Fock [1946] given as follows, 

wi(t) = i Vrr [Ait) - iBi(t)l (34) 

one obtains the following expression for the field on the surface 

1 f x_ir h{t) 
u(y) = ir- ~ ~^7: ex P Uks + iyt/d} dt, (35) 

ZTTl J -yc-jcW^t) 

where - t 

h{t) = 2i\/Trf3-\- w\{x)dx, (36) 

Jo 

and , m 2 2^ 3 /X/? 2 V /3 

d= ~ir = {-v) ' w 

The parameter d is the distance that Fock [1946] gives for an estimate of the width of the 
penumbra region produced by a plane wave incident upon a parabolic surface. For the case at 
hand the parameter d may be thought of as being an estimate of the width of the transition region. 
For y sufficiently large (y the order of d), the integral representation may be expressed in terms 
of the residues at the zeros of w[(t), giving the result expressing the field in terms of the creeping 
waves 

u{y) = e ik * 2 Ci exp [iytjd] (38) 

i=i 
where 

C/=[2iV^/3 + j l wi(x)dx]l[tiwi(ti)] ( 39) 

with ti being the roots of w[ (t) = 0. 

For computational purposes, it is best to use the following relationship involving the Airy 
integrals: 

wi(t) = e™>« 2\Tn Ai (fe +2 ^ 3 ). (40) 

Setting ti = e i7Tl3 p h (38) and (39) become 

u(y) = ^Ci exp {iks + iypi e^/d} (41) 

i=i 

with 

Q= [1/3+ f 'AA-x)dx]l[fiM-fid]. (42) 
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The values of /3* (the roots of A\{— /3) = 0) and A\{— f3i) can be obtained from tables given by Logan 

[1959]. The integrals Ai(—x)dx can then be interpolated from the National Bureau of Stand- 
Jo 

ards Tables [1958]. Proceeding in this manner one can easily compute the coefficients Ci. 

Particular values of the appropriate coefficient involved are given in table 1 for the first three roots 

(Z=l,2, and 3). 



Table 1. 


Roots and residue co- 
efficients 


/ 


0i 


A A- /3,) 


AA-x)dx 


Ci 


i 

2 
3 


1.0188 
3.2482 
4.8201 


0.5357 

-.4190 

.3804 


0.4757 
.7008 
.6511 


1.4823 

-0.7598 

.5369 



Away from the locally parabolic section the field on the surface given by (41) may be generalized 
to include a general convex surface of varying curvature by using the fact that Ci are the launching 
factors of the creeping waves. Thus using the general expression for the creeping waves [Keller 
and Levy, 1959] one has 



u(s)- 



k{s 



»l/6 



k(0) 



^Ci exp i f ks+ e^Pi QV /3 [V 3 ds\ 



(43) 



where k(s) is the curvature at a distance s along the surface measured from the z-axis. 

3. Convex Section in the Illuminated Region 

Having considered the case where the incident wave is propagating along the planar section 
with the parabolic section in the shadow region, the case where the plane wave is incident to 
the planar section with the parabolic section in the illuminated region, will be briefly considered. 
The same coordinate system and surface will be used as in the part 2, but the incident plane 
wave will have the form 

E = Le- i ^V/x /€o- 
Again (2) will be used to describe the total magnetic field on the surface, but with 

u (p) = e- ik v. 
The fields on the surface will be assumed to have the general form 

y > u(y) = e~ ik v I+(y) + e iks K+(y) 



y < u(y) = e~ ik y /_(y) 



(44) 



(45) 



where the second term in (44) is the wave reflected at the junction. Since the surface is smoothly 
joined, this reflected wave is of lower order magnitude and to the first approximation can be 
neglected. 

Using (3), but with the nonhomogeneous term exp(iA:y) replaced by 2 exp(— iky), together with 
(4) and (5) and relations (44) and (45), one obtains the following, holding for y< 0, 



I-iy) 



4/? Jo 



dtH[ l Hkr) 



(? 2 -2ty) 
r 
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{/+(*) €«*<»-<>+ K + (t)e^» +,) } 



(46) 



where r is prescribed by (9a). Using the relationships given by (15) and (16), together with 

I-(y)=J-(£) 

I + (y)=J+({), 

and taking the asymptotic limit as /n— »°°, (46) reduces to 



/_<0-2— J-yJ-J ^77^)^ eXP l8(r 



r-^)r + 



(t)[1+0 



(47) 



In a similar manner one can proceed and obtain the following integral equation holding 



for y > 0, 



/+tf)=2-- 



-in/ 4 



7T Jt 



*/>)(t-0^«pJJ(t+0»(t-0 



(48) 



The solution to this equation at a point £ depends only on the values of the field at points in the 
illuminated region such that r ^ £, hence is independent of the planar section. Thus, as is ex- 
pected, 7 + (£) is the solution that would be obtained if the complete surface were parabolic. Thus 
to a first approximation in kR, we have for y > from [Fock, 1946] 



Ufy 



v£L* 



exp [— Mj\2~ 



1/3 _ 



•jgVfl 



«\'(t) 



(49) 



where the contour Ti runs from infinity to the origin along the ray arg £ = 2/3 77, and from the origin 
to infinity along the positive real axis of t. 

Combining (47) and (49), /-(f) can be given explicity. The main attention will be focussed 
upon the asymptotic evaluation of /-(f) as |f| becomes large. This will then yield an asymptotic 
relationship for u(y) on the flat segment a long distance from the join. Since /+(r) has the following 
asymptotic expression for large r, 



/ + (t)~2+o 



(50) 



the integral in (47) will be broken up into two portions. The first portion will be of the form 



e -mi4 |2 f» (r 2 -2rf) 



8(r-a 



U4r)-2} 



which reduces in the limit |f | — » «> 

e ~i7r/4 / 2 \l/2 



2 \ i/2 r * 

=h) J. «*iw>-«- 



(511 



The second portion is of the form 



e ~ml4 J2 f: 



<(~2 



dr 



(T 2 ~2Tg) 



:8(t-0. 



i.e., it is obtained from the second term on the rightdiand side of (47) by replacing J+(t) by its 
asymptotic value 2 corresponding to the physical optics approximation for the current. The above 



774-625 0-65— 5 
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expression has the asymptotic value (obtained by setting r = z(— 8£) 1/4 ) in the limit as |f | — > oo : 

- e-™' 4 J- J exp [*z 4 ]dz 2 + 1 f | " 3 / 4 , (52 ) 

which reduces to — 1 + 0|£|~ 3/4 . On using (51) and (52), one obtains the following expression for 

y — >00; 



u(y) ~ e- ik «[l-(2TT | £ | j)" 1/2 \T[J + (r)-2]dT. 

Jo 



(53) 



The first term corresponds to the incident field. The second term representing a traveling wave, 
is the order of |A:y|~ 1/2 . 

4. Discussion 

The results in the preceding two sections can be used separately or together. For instance 
the surface could be comprised of three segments, two parabolic convex sections smoothly joined 
to and enclosing a flat section. In this case one convex section would be in the geometric optics 
illuminated region and the other in the shadow region. If the length L of the flat section is a suf- 
ficient number of wavelengths long, then the field distribution on the surface a long distance 
away from the first joint (the join of the illuminated convex section with the flat section) would 
be given by (53). Thus the field incident upon the second joint would have amplitude prescribed 
by (53) with the appropriate value for f (i.e., |y|=L). The field distribution in the shadow 
region would be prescribed by (35) or (43), modified by the appropriate factor associated with the 
value of the amplitude and phase of the incident field at the second joint. 

Even though the results were derived for cylindrical surfaces, they could be applied locally 
to a general convex surface provided that the radius of curvature transverse to the direction of 
incidence was very large compared to wavelength, and the shadow boundary coincides with the 
join. A special case would be for nose-on incidence to a body of revolution, formed by rotating 
the profile of the surface about an axis which is parallel to the flat portion of the profile. The 
results of the second section could be employed to obtain the expression for the creeping waves 
on the spherical base of a cone-sphere shape where the incident wave would be a wave launched 
from the tip of the cone. 

Finally a few comments should be made regarding the horizontally polarized case. In place 
of (2), the following is employed for the horizontal polarization: 

du duo 1 f dG du 70 



Lf. 



drip drip 2tt J drip driQ 

where u and — are proportional respectively to the electric intensity and the surface current (both 
on 

vectors directed in the z direction). For the geometry treated first wherein the incident wave is 

a plane wave traveling along the flat portion towards the join with the convex section in the shadow 

zone, the incident field on the surface I associated with the term 2 — — ) is zero. Hence the field or 

\ dn/ 

current is zero everywhere. A nonvanishing current distribution will be obtained only if the 

source is included (such a source could be that produced by an edge). For the other geometry 

where the convex section is in the illuminated zone, the analysis is similar as that for the vertically 

polarized case. However a long distance away from the junction on the flat portion, the current 

distribution is just comprised of a traveling wave which decreases in amplitude. 
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5. Appendix A 

The following integral 

F(p) = j* exp {-p*-f } f" (I)'' 2 exp [i (|) 3 r} drdz (A .l) 

will be evaluated in terms of Airy integrals. First differentiate the above expression for F(p) twice 
with respect to p and integrate by parts giving 

- l -F"(p)-pF(p) = -^ exp (iT*)dT + ^^ exp ("P*-|Q * 1/2 ^ 

where the prime indicates differentiation with respect to p. The first term on the right-hand side 
is the Fresnel integral and is easily evaluated, and the second integral can be expressed in terms 
of Airy integrals [Weston, 1962]. The right-hand side thus becomes 

- Vf e'"< + 1 »*»«** j- q {AiqUiq) - iB t (q)]} (A .2) 

where the argument q of the Airy integrals are 

q = i2 l > 3 p. (A.3) 

Set p = t 2-m e -ini6 (A4) 

F(p)=F(t). (A.5) 



The resulting differential equation for F(t) becomes 

<PF(t) 



dt 2 
^here 



tF(t)=f(t) (A.6) 



fit) = V^r 2-WV" 8 " 2 [- 1 + 3m ^ {AM.AM + «#)]}] (A.7) 

is derived upon using the relations 



^9)=-|e*^ 3 [^)+ifi^)]l 



^ i (9)-Jfi.(?) = -2e 2 ' r '/^ i (r) 



(A.8) 



The solution of the homogeneous portion of the differential equation for F(t) are the Airy 
integrals; hence the complete solution is given by 

F(t) = dAit) + C 2 Bi(t) - 7T JJ [A i (t)B i (x)-Ai(x)B i (t)V(x)dx (A .9) 

where C\ and C<l are constants yet to be prescribed. The integral can be reduced using the fol- 
lowing relations 
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+ A*B' t \dx 



3 \* A\B\dx = [ t §2A i A' i Bi + -A i 
= -[ t A i (x)dx + [A*B i y 
3 f A\AiBidx = -- [* Ai{x)dx + [A*Bi]< 

JO 77 Jo ° 

3 {' B'.B i A i dx = - P Bi(x)dx + [B*Ai]< 

JO 7Tjo ' ° 



giving 



#(t) - 2 1 / 3 7r 3 / 2 e--/ 12 [^ / (f) -f £B^)] f X ^(*)<fe + GA(*) + C 2 Bi(t) (AJ0 ) 

where Ci and Cz are two new constants. 

The next problem is to determine these constants. We shall first consider the behavior as 
| ^| — > oo with |arg t\ < 77-/3. In this case 

m ~zfa> eip {-H (A ' n) 

m ~V^f. exp (f '"1 < A - 12 > 

and F(t) has the asymptotic value 

^)~Se« + ^exp{|^}. (A . 13) 

The term with coefficient C\ is subdominant and can be ignored. Returning to the original expres- 
sion for F(p) we have, on integrating by parts, 

F{p) ~ij\ir> dT+0 (±y (A . 14) 

Thus F(t) must have the form for large \t\ 

^-f^^' 12 (A.15) 



and it then is seen that the coefficient C2 must vanish. 

77 77 

To evaluate Ci, we will again let | f | — > 00 but in the sector — < arg t < 277/3 (i.e., — < arg 

p<f). In this case we shall set 

I Ai(x)dx= I Ai(x)dx— Ai(x)dx. (A. 16) 
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Setting t = re 2nil3 with |arg r| < fir/3, 

the asymptotic behavior of ^i(0 and Z?i(0 becomes 



Ai{t) = -\e^[MT) - iBi(r)] ~ * exp {tt*/3 + 2/3T 3 ' 2 } 

Z Z \/7r T ' 

^-(;) + iB 1 (t) = -2e 4iri ' 3 A i (T) ~ 7= , „ exp {777/3 - 2 /3r 3 ' 2 }. 

yirr* 

It can be shown that 

/„^ lWA ~2^^ eXP {f T3/2 }' ■ < A - 17 > 

hence we have (on neglecting the subdominant terms) 

FU) ~^^ + ^& exp <™/3 + 2/3t*} (A . 18) 

and thus it is seen that C\ must vanish also. 
We have for the final result 



F(t) = 2 1 l B v 3 l 2 e- iv l 12 [A i (t) + iBi(^ \' A&x)dx 

The definite integral may be expressed in the form 

foe Q-l/2 foe 

Ai(x)dx = K m (y)dy 

JO 7T Jo 

involving the modified Hankel function, which reduces to the result 

f« 1 

Ai(x)dx = -- 
Jo o 



(A.19) 



(A.20) 



For convenience, (A.19) will be expressed in terms of Airy integrals with argument q given by 
(A.3). Thus from (A.3) and (A.4) we have 

t = e 2 ™l 3 q 

and using relationships (A. 8), (A.19) becomes 

F(p) = &W****AM) {!+! | f \Ai(x) - iBi(x)]dx} , (A.21) 

with q = -ip2 x l\ 
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